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Abstract 

The relaxation kinetics of cavity polaritons by scattering with thermal acoustic phonons is studied 
within the rate equation approximation. Numerical results show that a suppression of the bottleneck 
£S| . of lower polariton states occurs at high polariton densities. We have found that the long decay time 

of the photon-like polaritons, the thin width of the embedded quantum wells and the small value of 
' exciton-cavity detuning are favorable for the suppression of the bottleneck. 
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CT> ! 1 Introduction 

Semiconductor microcavities (MC) with embedded quantum wells have found a growing interest in recent 
years. The eigenstates of the systems as a result of the strong coupling between the MC photon and 
quantum well exciton at the same in-plane wave vector are called cavity polaritons. The Rabi splitting 
and the dispersion of these cavity polaritons are readily observable using reflection, absorption, photo- 
\ luminescence measurements [1-3]. Recently, many new phenomena have been intensely studied in MC: 

f***- ■ super-linear behaviour in emission intensity [4-6] , stimulated polariton amplifier [7-8] and suppression of 

the relaxation bottleneck 0. To the best of our knowledge, Tartakovski et al. were the only group, 
which has clearly observed this suppression of the relaxation bottleneck of lower polariton states at higher 
powers. 

In this work, we study the relaxation kinetics of cavity polaritons by scattering with thermal acoustic 
phonons. Here we will present numerical solutions of the rate equations and show that a suppression of 
I ■ the bottleneck occurs at sufficiently high concentration of the polaritons. 
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H ; 2 Relaxation kinetics of cavity polaritons 



Within the two-coupled band model and considering only the Is state of the heavy-hole quantum well 
exciton, the energy for the upper and lower cavity polariton branches is given by the usual solutions of 



the polariton dispersion equations as given in Rcf. [|10| 
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where 
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El = E% + — (3) 

k 2(m e +m h ) v ; 

Here Eft is the fundamental exciton energy, m e (ma) the electron (hole) in-plane mass, n cav the cavity 
refraction index, f2 the Rabi splitting. Fig. 1 shows the cavity polariton dispersion of GaAs quantum 
wells embedded in cavity for a detuning 6 — Eq — Efi = —2 meV, a Rabi splitting nil — 4.8 meV, and 

Hc.ai) — \/ ^r.n.u 3.07. 
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In the following we will use the semi-classical Boltzmann kinetics to investigate the cavity polariton 
relaxation by acoustic phonon scattering. The rate equation for the distribution ffc of branch i with the 
in-plane wave vector k and energy E\ of the cavity polariton can be written as [10-11] 
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~ is the radiative recombination rate and G\.(t) = A l k H{t) the generation rate. The transition rates 

Will — >/ due to deformation potential scattering are weighted with exciton Hopfield coefficients x\ and 
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Here L z is quantum well width, V the sample volume, u the longitudinal sound velocity , p the mass 
density of the solid, D e and Dh the deformation potentials. A(x) and G(x) are the Fourier transform of 
Xe( z e) and F 2 (p), respectively. Note that the Is exciton wave function has been taken in the simplest 
form 
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4>x S ( r e,r h ) = F(p)xe{ze)Xh{z h )—=e 



(11) 



Expression (||) contains transition rate due to the absorption of phonon oc iVg for E > 0, and due 
to emission of phonon oc (N^ El + 1) for E < 0, where N? >Q = (e pE - l) 1 is the thermal phonon 
distribution function. 

In the following wc will present the results of the cavity polariton kinetics on the lower branch assuming 
an isotropic distribution of polaritons in k space. For the decay rate ~, we use the same approximation 



(4T 



for < k < k c 



6 x 10 4 cm 1 : 



given by Bloch and Marzin in GaAs quantum well |ilf| : i 
~ = 4j- for k cav < k < k ra d — n " a t " E ° = 2.1 x 10 5 cm -1 ; and ~ = for k > k ra d, (r c and t x are the 
radiative life of photon and exciton, respectively, ci the photon Hopfield coefficient). 



3 Numerical results for kinetics of cavity polaritons 

In this section, rate equations are numerically integrated using the following parameters for MC with 
embedded GaAs quantum wells: [12-13] _Eq = 1.3955 eV, m e = 0.067mo, mh = 0.13mo, hfl = 
4.8 meV, a b = 7.5 nm, D h = 12 eV, D e = -7 eV @, u = 4.81 x 10 3 m/s, p = 5.3 x 10 3 kg/m 3 , n cav = 
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3.07, t x — 200 ps. Concerning the generation rate we choose for A\ a Gausian distribution of exciton 
on the lower branch with a width AE = 0.1 meV centered at a momentum fco = 5 x 10 5 cmT x and 
H{t) = Hq th (y) with T — 20 ps. Unless otherwise stated we will choose L z = 6 nm, T p h on = 2 K, 
5 = E§-E$ = -lmeV. 

In Fig. 2 we show the distribution function fk(t) of cavity polaritons on the lower branch for a rather 
long radiative photon decay time t c = 50 ps at n — 7.8 x 10 9 cm" 2 . We clearly see the presence of 
the bottleneck effect similar to the bulk case [15-16]. However, if we increase the pump intensity, a 
suppression of the bottleneck occurs as can be clearly seen in Fig. 3 at n = 5.1 x 10 10 cm" 2 . These 
findings agree rather well with recent measurements performed in GaAs cavity by Tartakovski et al. || 
using nonresonant cw laser excitation. They found that the bottleneck has been strongly suppressed at 
higher power laser excitation. This feature is specific for cavity polaritons. We have also carried out 
the calculations for the ortho-exciton polaritons in bulk Cu^O |l6| and have found that the bottleneck 
persists up to very high density around 10 18 cm -3 . The main difference in the dispersion relation between 
cavity polaritons and bulk polaritons is the slope of E^ versus k in the bottleneck k-region. The very 
high value of the slope of the bulk polaritons (70 times higher than the slope of the cavity polaritons) 
is very unfavorable for the suppression of the bottleneck even at very high densities. On the other 
hand, one can lower the slope in the cavity polaritons by decreasing the absolute value of the negative 
detuning 6. From expression (JoJ) of A(x) we can see that the scattering rate is roughly proportional 
to (sink^/Lz} . Therefore, quantum wells with thin width will be favorable for the suppression of the 
bottleneck. For L z = 3 nm, 8 — —0.25 meV and t c = 50 ps the bottleneck starts to be suppressed 
at densities, which are roughly 3 time lower than those for L z — 6 nm, 8 = — 1 meV and r c = 50 ps 
. For a more realistic r c one needs a higher density to compensate the fast decay of the photon-like 
polaritons. Recently, Senellart and Bloch have given the photon decay time r c a value of 6.6 ps |l7|. In 
Fig. 4 we display the integrated distribution function ( dtfk{t)) versus the wave number k for different 
polariton densities for r c = 6.6 ps, L z — 3 nm, 8 = —0.25 meV at 2K (Fig. 4a) and 15K (Fig. 4b). These 
results show a clear process of the relaxation kinetics: a nonlinear increase of low k cavity polaritons with 
increasing polariton densities (or excitation laser powers) leads to the suppression of the bottleneck. The 
suppression is slightly more favorable at T=15 K. For density n > 10 10 cm~ 2 , one should take into 
account the exciton-exciton scattering and also go beyond the rate equation approximation. However, we 
expect our results remain qualitatively valid even in a better treatment of polariton relaxation kinetics. 

In conclusion, we have presented detailed cavity polariton relaxation kinetics of lower branch due to 
deformation potential scattering by accoustic phonons. In agreement with recent experimental results by 
Tartakovski et al. |J we found that suppression of the bottleneck is found at high polariton densities. 
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Figure captions: 

Fig. 1: Dispersion of cavity polaritons with embedded GaAs quantum wells for S = — 2 meV, TiVL = 
4.8 meV. 

Fig. 2: Resulting distributions of cavity polaritons on lower branch versus time for 6 = — 1 meV, L z = 
6 ran, r c = 50 ps, n — 7.8 x 10 9 cm -2 at a bath temperature of 2 K. 

Fig. 3: Suppression of the bottleneck at n — 5.1 x 10 10 cm -2 for 6 = — 1 meV, = 6 nm, t c = 50 ps 
and T = 2 if . The inset shows total polariton density. 

Fig. 4: Resulting integrated distribution function versus wave number k for different polariton den- 
sities at 2K (a) and 15K (b) for S = —0.25 meV, L z = 3 nm, t c — 6.6 ps. 
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Fig. 1: Dispersion of cavity polaritons with embedded GaAs quantum wells for 5 = 
-2 meV, htt = 4.8 meV. 




Fig. 2: Resulting distributions of cavity polaritons on lower branch versus time for 
S = — 1 meV, L z = 6 nm, r c = 50 ps, n = 7.8 x 10 9 cmT 2 at a bath temperature of 2 K. 




Fig. 3: Suppression of the bottleneck at n — 5.1 x 10 10 cm 2 for 5 = — 1 meV, L z = 
6 nm, r c = 50 ps and T = 2 K. The inset shows total polariton density. 
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Fig. 4: 

Resulting integrated distribution function versus wave number k for different polariton 
densities at 2K (a) and 15K (b) for 5 = —0.25 meV, L z = 3 nm, t c = 6.6 ps. 



